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ABSTRACT 

SONYC - Substellar Objects in Nearby Young Clusters - is a survey program to investigate the 
frequency and properties of substellar objects with masses down to a few times that of Jupiter in 
nearby star-forming regions. Here we present the first results from SONYC observations of NGC 1333, 
a ~ 1 Myr old cluster in the Perseus star-forming complex. We have carried out extremely deep optical 
and near-infrared imaging in four bands (i', z', J, K) using Suprime-Cam and MOIRCS instruments 
at the Subaru telescope. The survey covers 0.25sqdeg and reaches completeness limits of 24.7 mag 
in the i'-band and 20.8 mag in the J-band. We select 196 candidates with colors as expected for 
young, very low-mass objects. Follow-up multi-object spectroscopy with MOIRCS is presented for 53 
objects. We confirm 19 objects as likely brown dwarfs in NGC 1333, seven of them previously known. 
Nine additional objects are classified as possible stellar cluster members, likely with early to mid M 
spectral types. The confirmed objects are strongly clustered around the peak in the gas distribution 
and the core of the cluster of known stellar members. For 11 of them, wc confirm the presence of disks 
based on Spitzer/IRAC photometry. The effective temperatures for the brown dwarf sample range 
from 2500K to 3000K, which translates to masses of ~ 0.015 to 0.1 M Q , based on model evolutionary 
tracks. For comparison, the completeness limit of our survey translates to mass limits of 0.004 M Q 
for Ay < 5 mag or 0.008 Mq for Ay < 10 mag. Compared with other star-forming regions, NGC 1333 
shows an overabundance of brown dwarfs relative to low-mass stars, by a factor of 2-5. On the other 
hand, NGC 1333 has a deficit of planetary-mass objects: Based on the surveys in a Orionis, the Orion 
Nebula Cluster and Chamaeleon I, the expected number of planetary-mass objects in NGC 1333 is 
8-10, but we find none. It is plausible that our survey has detected the minimum mass limit for 
star formation in this particular cluster, at around 0.012-0.02 Mq. If confirmed, our findings point to 
significant regional/environmental differences in the number of brown dwarfs and the minimum mass 
of the Initial Mass Function. 

Subject headings: stars: circumstellar matter, formation, low-mass, brown dwarfs - planetary systems 



1. INTRODUCTION 

The origin of the stellar Initial Mass Function (IMF) 
is one of the major issues in astrophysics. The low-mass 
end of the IMF, in particular, has been the subject of 
numerous obser vational and theoret ical studies over the 
past decade (see lBonnell et al.ll2007|) . While higher mass 
stars exhibit the well-known Salpete r shaped mass func- 
tion (dN oc m~ Q dm, with a = 2.35. lSalpeterlll955D . the 
power-law slope becomes significantly sh allower below 
- 0.5 M Q (IKroupal[200lt IChabrierl[200l . In the sub- 
stellar regime, t he coefficient a is mostly found to be 
around 0.6 (e.g. lMoraux~et al. 2003; Lopez Mart i et al.l 
l200i iMoraux et al.l[2007t ICaballero et al.ll200l . While 
the higher-mass domain is thou ght to be mostly 
forme d through fragmentation (e.g. lPadoan fc Nordlundl 
I2002D and /or accret ion onto the protostellar core 
(|Bonnell fc Batej|2006l e.g.), in the low-mass and substel- 
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lar regime additional physics is likely t o play an impor- 
tant r ole, e .g. dynamical i nteractions dBate et al.l 120031 : 
iBatd I2009D . turbulence (|Padoan fc Nordlundl |200¥ . 
photoerosion in the radiatio n field of bright stars 
(|Whitworth fe Zinneckerl 120041 ). or tidal shear in the 
gravi tational potential of a stellar cluster (|Bonnell et al.l 
l2008h . 

To continue collapsing, an object or region has to cool 
as quickly as it is being heated by the conversion of grav- 
itational into thermal energy. As radiative cooling is re- 
lated to the opacity, this condition defines the opacity 
limit of fragmentation, which translates into a minimum 
mass for star-like sources, predicted to be between 0.001 
and 0.01 Mq (e.g. lLow fc Lvnden-Belll 119761: iReesI 119761 : 
lBossll2001b iBatel 120051: IWhitworth fc Stamatellosll2006h . 
Thus, isolated substellar objects in star-forming regions 
can have planetary-like masses. Furthermore, current 
numerical simulations of star formation build a signifi- 
cant fraction of brown dwarfs (BDs) by fragmentation 
in pr otostellar disks and subsequent ejection from th e 
disk (jBate et al.1 120031 : IStamatellos fc Whitworthll2009h . 
Thus, a fraction of the substellar population may have a 
similar origin to massive planets. Hence, using the term 
'iso lated planetary-mass object' (IPMOs or planemos, 
see lBasri fc Brownll2006l ) for objects with masses below 
the Deuterium burning limit (~ 0.015 M©) makes some 
physical sense, as this is the mass regime where star and 
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planet formation overlap. 

The frequency and physical properties of the ob- 
jects at the extreme low-mass end of the spectrum are 
oorly known . Following the discovery of b rown dwarfs 
Rebolo et all 119951 : iNakajima etaD Il995l ). surveys of 
star-forming regions, young clusters, and the field have 
revealed a rich population of substellar objects, trigger- 
ing a wealth of follow-up studies to characterize their 
physics (see B asrH 120001 ). Populations of planetary- 
mass objects have been i dentified with masses down 
to 0.005 - 0.01 My (e.g. IZapatero Osorio et alj 120001 : 
iLucas fc Rochd l200Q[ ) . So far, there is no evidence for 
a cut-off in the mass function, or in other words, the 
opacity limit has not been observed yet. 

As of today, however, only two regions have been sys- 
tematically surveyed down to masses significantly below 
the Deuterium burning limit: erOrionis cluster and the 
Orion Nebula Cluster (ONC). At a distance of ~ 400 pc, 
any follow-up for planemos in Orion beyond photom- 
etry and very low-resolution spectroscopy is challeng- 
ing. A few more objects in the same mass regime have 
been identified based on mid- infrared detections from 
ISO or Spitzer (TTesti et all l2002t ILuhman et all 120051: 
I Alters et al.l l2006t iLuhman fc Muenchl l2008fh This ap- 
proach, however, will by definition only find sources with 
circumstellar material and thus yield incomplete results. 
The unbiased census for the members of the Chamaeleon 
I (Chal) star-f orming region is currently complete down 
to - 0.01 M Q (|Luhmanll2007h . and has yielded a handful 
of planetary-mass members. 

SONYC - Substellar Objects in Nearby Young Clus- 
ters - is a new, unbiased attempt to establish the fre- 
quency of low-mass brown dwarfs and planetary-mass 
objects as a function of cluster environment and to pro- 
vide the groundwork for detailed characterization of their 
physical properties (disks, binarity, atmospheres, accre- 
tion, activity). The project is based on extremely deep 
optical and near-infrared imaging and follow-up spec- 
troscopy using 8-m class telescopes, aiming to detect the 
photospheres of the objects, which is a prerequisite for 
an unbiased survey. The observations are designed to 
reach limiting masses of 0.003 Mq, and thus are meant 
to probe the opacity limit of star formation (see above). 
In each region, we aim to cover at least ~ lOOOarcmin 2 . 

In the first paper for this project, we report on the 
survey in the cluster NGC 1333, part of the Perseus 
star-forming complex, carried out with Suprime-Cam 
and MOIRCS at the Subaru tel escope. With an age 
of ~ lMyr, a dista nce of 300 pc (|de Zeeuw et alJH999t 
iBelikov et al.l [2002). and moderate extinction, it is fea- 
sible to reach completeness limits of I = 25 mag and 
J = 21 mag, corresponding to masses for cluster mem- 
bers of about 0.0 03 Mp,, according to the COND03 evo- 
lutionary tracks (jBaraffe et al.l l2003). Furthermore, the 
cluster is compact and can thus be covered efficiently. 

Near-infrared survey s in NGC 1333 have been carried 
out, a mong others, by lStromet al.1 (|1976h . lAspin et al.1 
(|1994h . and lLada etatl (|1996t ). establishing the presence 
of an embedded stellar cluster and a site of ongoing 
star formation with more than a hundred YSO s down 
to K = 14.5 mag. More recently, IWilking et~afl (|2004l ) 
have identified at least nine substellar objects based on 
near-infrared imaging over 133arcmin 2 with a limiting 
magnitude of K = 16 mag and follow-up spectroscopy. 



They find a disk frequency of 74% and an upper limit 
to the mass function exponent of a < 1.6 in the sub- 
stellar regime, consis tent with prev i ous s tudies. Build- 
ing upon this work, iGreissl et al.l |2007) have carried 
out deeper HST/NICMOS observations of a limited area 
(3 arcmin 2 ), probing for the first time the planetary-mass 
regime. The near-infrared survey by lOasa et"aT1 (2008) 
is similarly deep (K — 18 mag) and covers 25 arcmin 2 , 
but does not include follow-up spectroscopy. In ad- 
dition, NGC 1333 has been covered by numerous X- 
ray, mid-infrared, sub-millimeter, and millimeter sur- 
veys, for example in the framework of the Spitzer 'Cores 
to Disks' Legacy Program, resulting in detailed stud- 
ies of the cloud and cluster stru cture as well as ac - 
cretion and outflow activity (e.g. iGetman et all | 200 
| Preibischll2003t lEnoch et alJl2006b I Jorgensen et aljboO 
Wal sh et alj|2007t iGutermuth et al.ll2008h . 

This paper is structured as follows: We present our 
observations and data reduction procedures for photom- 
etry and spectroscopy in Sect. [2] In Sect. [3] we discuss 
the photometric selection and spectroscopic verification 
of very low-mass objects in NGC 1333. The properties of 
our sources, including evidence for youth, are analysed in 
Sect, m The low- mass end of the IMF based on our new 
findings is discussed in Sect. OH We give our conclusions 
in Sect. El 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Optical imaging and photometry 

We observed NGC 1333 in the SDSS i' and z' filters 
using the Subaru Prime Focus Camera (Suprime-Cam) 
wide field imager (jMivazaki et alJl2002h on 2006 Novem- 
ber 17. The conditions were photometric with typical 
seeing ranging from 0.45-0.63". Suprime-Cam is a mo- 
saic camera which utilizes 10 CCDs arranged in a 5 by 
2 pattern giving a total field of view of 34' x 27'. As 
NGC 1333 is quite compact, the entire cluster was ob- 
served in a single pointing. A ten point dither pattern 
was carried out to eliminate gaps between the CCDs and 
to correct for bad pixels. We observed this pattern in 
each filter six times, where the individual images have 
exposure times of 60 sec, for a total integration time of 
3600 sec in each band. 

A standard reduction was performed for each individ- 
ual CCD chip. Routines from the S uprime-Cam data 
reduction softwar e package SDFRED (jYagi et al.l 120021 : 
lOuchi et al. 2004) were used for overscan subtraction and 
flat fielding. In both filters, the images show significant 
extended structure from the molecular cloud emission. 
Sky subtraction was performed by fitting a Gaussian to 
the peak of a histogram of pixel brightnesses and sub- 
tracting the centroid. Additional routines from SDFRED 
were used for distortion correction, bad pixel masking 
and image combination. Only one of the ten CCD chips 
(w67cl) showed noticeable fringing. To reduce each chip 
consistently, and because the 10 point dither pattern 
mostly removes the fringes, we did not carry out a fringe 
correction. The world coordinate system was calibrated 
against sources fro m the 2MASS point source catalog 
(jCutri et al.l |2003j) with J > 13 using the msctpeak pro- 
gram from the IRAF package MSCRED 5 . The typical 
fitting residuals were of order 0.25". 

5 IRAF is distributed by the National Optical Astronomy Ob- 
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Fig. 1. — i'-band histogram of the objects in our photometric 
catalogue with i'- and z'-band data. The peak at 24.7 mag indi- 
cates the completeness limit of the survey. The faintest objects in 
the survey are found at i > 26 mag. The dashed and dash-dotted 
lines show the histogram for equally sized samples of blue and red 
objects, illustrating that the completeness limit does not signifi- 
cantly change with color, although we are slightly more sensitive 
to red objects. 

Sources were identified from each CCD chip us- 
ing the Source Extractor (SExtractor) software package 
(|Bertin fc Arnouti 119961 ). For object extraction we re- 
quired at least 5 pixels to be above the 3a detection 
limit. SExtractor automatic aperture fitting photometry 
routines were used to calculate the flux of each source. 
We rejected objects within a few pixels of the edge of 
the image, elongated objects (a/b > 1.2) and objects 
that were not present in both the i' and z' source lists. 
The rejection criteria were chosen conservatively, to min- 
imize the number of spurious detections in the photom- 
etry database. We visually checked the detection algo- 
rithm and confirmed that the number of point-sources 
that are missed as well as the number of spurious sources 
that are included is minimal. The final optical catalogue 
has 7757 objects. 

We measured the chip-to-chip zeropoint offsets from 
the median fluxes of domeflat images. The absolute ze- 
ropoint for the mosaic was derived from observations of 
the SA95 s tandard field, whic h contains SDSS secondary 
standards ( Smith et al.l 2002). To avoid saturation, the 
images of standard fields were defocused. The fluxes of 
five standard stars were extracted using aperture pho- 
tometry with 50-pixel radius. We derived a mean ze- 
ropoint of 27.80 ± 0.01 in F- and 27.12 ± 0.02 in z'- 
ba nd. These valu e s agr ee with the zeropoint measured 
bv lMivazaki et ail (|2002l ) within 0.12 magnitude in i' and 
0.07 in z'. The airmass difference between standard fields 
and science fields was ±0.1, i.e. the differential extinction 
is negligible. 

From the histogram of the magnitudes (see Fig. [1] 
for the i'-band data) we derived completeness limits of 
our survey: In the i'-band, we are complete down to 
24.7 ±0.2 mag, in z'-band down to 23.8 ±0.2 mag. These 
limits are robust against changes in color, as demon- 
strated in Fig. [T] The limiting magnitudes of our sur- 
vey are about 26.5 mag in i' and 25.5 mag in z', which 

servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



are close to the sensitivity limits det ermined during 
the S uprime-Cam commissioning phase (Mivazak i et al.l 
l2002h . 

2.2. Near-infrared imaging and photometry 

NGC 1333 was observed with the Multi-Object In- 
fraRed Camera and Spectrograph (MOIRCS) mounted 
on the Subaru Telescope, in the MKO broad-band fil- 
ters J and Ks, during the nights of 2008 December 4-6. 
MOIRCS uses two CCDs providing a total field of view 
of 4' x 7' with a spatial resolution of 0.117 arcsec/pixel. 
NGC 1333 was covered in 24 pointings in both filters, 
using a 6 or 8-point dither pattern. In the J-band, the 
typical exposure time per pixel was 600 sec, in the ri- 
band 300 sec. In addition, the center of NGC 1333 was 
covered in a smaller series of 13 sec short exposures, to 
obtain photometry for the brightest sources that saturate 
in the deeper exposures. 

Data reduction was performed using a modified version 
of the SIMPLE-MOIRCS package 6 , which is written in 
IDL 7 and uses SExtractor. Background subtraction, flat- 
field correction and distortion correction were performed 
on individual images, before the images were co-added. 
The world coordinate system was calibrated against the 
2MASS point-source catalogue. Initial source identifica- 
tion is performed using SExtractor, requiring at least 3 
pixels to be above the 3a detection limit. Source rejec- 
tion was performed based on saturation and elongation. 
After running the identification and rejection routines, 
the sources were visually verified for each field. Photom- 
etry was performed using a fixed-aperture extraction rou- 
tine in IDL. To account for variable seeing between fields, 
an aperture correction factor is determined for each in- 
dividual field. To correct for difference in sensitivity be- 
tween detector 1 and 2, a flux offset factor is determined 
from the average level of the sky flux in both images, and 
applied to detector 2 for each individual field. 

The final near-infrared catalogue has 2360 sources. 
Calibration of the absolute zeropoint magnitude is per- 
formed by comparing the full catalogue of extracted 
sources to the 2MASS point source catalogue. The mean 
zeropoint is derived to be 25.08 ± 0.03 in J-band and 24.3 
± 0.04 mag in K-band. Completeness limits in each band 
are derived from the number distribution of the magni- 
tudes (see Sect. 12. ip . The imaging is complete down to 
18.0 mag in K-band and 20.8 mag in J-band. 

2.3. Multi-object spectroscopy 

We used MOIRCS again to carry out multi-object spec- 
troscopy (MOS) for 53 sources in NGC 1333, 36 of them 
selected based on the photometry in the four bands i', z', 
J, K (see Sect. [3]). The targets are covered in six multi- 
slit masks; the number of objects per mask was 5-13. 
Pre-imaging for the masks in the K-band was obtained 
in September 2008. We used slits which are 0.9" wide 
and 9-12" long, except for a few longer slits, used to ob- 
serve standard stars. The spectroscopy run was carried 
out in three nights on Dec 4-6, 2008. For each mask 
we integrated in total 50-60 min with the grism HK500, 
split in shorter exposures of 5 or 10 min, depending on 
conditions. Between the single exposures, the mask was 

6 http://www.aoc.nrao.edu/~whwang/idl/SIMPLE/MOIRCS/Doc/ 

7 Interactive Data Language 
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moved so that the targets shifted by 2.5" along the slit 
(nodding), to facilitate sky subtraction. The seeing dur- 
ing the science integrations was stable at 0.5-0.7". Be- 
fore and after the NGC 1333 masks, we observed AO 
stars through one of the science masks for flux calibra- 
tion and extinction correction; these exposures required 
de-focusing to prevent saturation. For each mask, we 
obtained series of domeflats with lamp on and off for 
calibration purposes. 

The MOS data was reduced following standard recipes 
for near-infrared spectra. In a first step, we subtracted 
the pairs of nodded exposures, which removes the sky 
background and the detector bias. These images were 
divided by a normalized Hatfield, which is the difference 
between the averaged lamp-on and lamp-off domeflats. 
Frames from the same nodding position were co-added. 
This gives us two final images for each chip and mask. 
In these frames, the signal-to-noise ratio for the spectra 
ranges from 20 for the faintest to 1000 for the brightest 
sources. These values have been measured in the central 
parts of the H- and K-band, in regions without strong 
skyline residuals. 

For the extraction of the spectra we used the apall rou- 
tine within IRAF. All objects from our photometric can- 
didate list are well-detected. For wavelength calibration 
we additionally extracted the spectra from the unreduced 
frames, which contain plenty of telluric OH lines. As the 
wavelength solution depends slightly on the position of 
the slit on the chip, the fit was done for each object sep- 
arately. For the HK grism, a second order polynomial fit 
covering about 20 skylines gave a typical rms of 2-3 A, 
well below the resolution. All spectra are smoothed with 
a ±25 A boxcar average and binned to 40 A per pixel. 
The two nodded spectra for the same object were coad- 
ded. The standard star spectra were treated in exactly 
the same way as the science frames. 

Near-infrared spectra show strong effects of the atmo- 
sphere - broad absorption features, mostly between the 
J/H and H/K bands, as well as sharp OH emission lines. 
In our data, subtracting the nodding exposures reduced 
the strength of the OH lines by 95% or better. The resid- 
uals are removed in the background subtraction which 
is part of the extraction. The corrections for telluric 
extinction and instrument response are usually done in 
one step, based on observations of a standard star with 
known spectral energy distribution. Our AO standard 
star spe ctra were div ided by a library spectrum of an 
AO star (|Pickleslll998h . obtained from the ESO website 8 . 
To remove residuals of the Hydrogen absorption features 
seen in AO spectra, we linearly interpolated from 1.54 
to 1.78 fim and from 2.09 to 2.34 /im, regions that are 
mostly free from additional telluric absorption bands. 

All science spectra were divided by one of these modi- 
fied standard spectra observed at similar airmass. Typ- 
ically, the airmass difference is < 0.1. Still, we find that 
the correction of atmospheric correction is not optimum. 
We attribute this mainly to the duration of our science 
exposures: While the standards were observed within a 
few minutes, the science exposures cover 1-2 hours. Vari- 
able sky conditions, as they have been present during the 
observations, hamper a good extinction correction. For 

8 http: / / www.eso.org/ sci / facilities / paranal/instruments/isaac/tools / 
lib/index, html 



our analysis we therefore focus on the spectral regions 
least affected by telluric absorption, i.e. 1.5-1.75 /im in 
the H-band and 2.1-2.3 /ito in the K-band. 

3. SELECTION OF NEW SUBSTELLAR MEMBERS IN 
NGC 1333 

3.1. Optical color-magnitude diagram 

The goal of the project is to identify new very low- 
mass members of NGC 1333. Due to their pre-main 
sequence status, plus extinction, these sources are ex- 
pected to occupy a distinct area in optical/near- infrared 
color-magnitude diagrams, on the red side of the broad 
cumulation of background main sequence stars. In Fig. 
[2] we show the (i', i'— z') diagram, which we used to 
identify the initial list of candidates. The cumulation 
of non-members around i' — z' = 1.0 mag is clearly 
seen. Previously identified very low-mass members of 
NGC 1 333 with spect r oscop ic co nfirmation from the pa- 
pers bv lWilking et alj (|2004h and lGreissl et all (|2007l ) are 
marked. 

Based on the locus of known members, we construct 
a selection box for this diagram, which encompasses the 
color-magnitude range where we expect cluster members 
with lower masses to be found (dash-dotted lines). The 
position of the box is chosen based on the colors of known 
substellar cluster members in the diagram. Apart from 
this empirical constraint, we aimed to select candidates 
with a broad range in colors, to maximize our chances of 
covering the cluster population. In particular, the choice 
of the selection box is independent from models. At the 
bright end, the box is limited by the onset of saturation 
at i' ~ 17.3 mag. The box extends down to our limiting 
magnitude of i' ~ 26.5 mag. The selection box contains 
196 objects with reliable photometry in i'- and z'-bands. 
From this initial sample, we select a representative list 
of 36 sources for follow-up spectroscopy, marked with +- 
signs in Fig. lH This sampl e inclu des 19 objects from the 
sample from lWilking et al. (2004), six of them previously 
spectroscopically confirmed as members. 

In the selection of candidates for spectroscopy, we ig- 
nore the bright end of the selection box i' < 18.5 mag, 
which contains about one third of the candidates. Clus- 
ter members with such magnitudes are likely to have 
masses around or above the substellar limit, and are 
therefore of less interest to us. Apart from that, we 
aim to sample the full range of colors and magnitudes 
in the box, to account for the uncertainties in estimat- 
ing the location of the cluster members in this diagram. 
Thus, the 'spectroscopic' sample does not show any bias 
in terms of colors or magnitudes in comparison with the 
full photometric candidate sample. In addition, there 
is no bias in spatial distribution: the full photometric 
sample shows a clear clustering around the cloud cores 
in NGC 1333, thus our spectroscopy fields are also dis- 
tributed in this area (see Sect I4.4|) . Since the density of 
young sources is higher in the regions close to the cloud 
cores, the contamination with field objects is expected to 
be larger further away from the cores. Thus, focusing on 
the cores is reasonable for minimizing contamination as 
well. 

3.2. Near-infrared color-magnitude diagrams 

We used the near-infrared photometry to verify the 
initial candidate selection. In particular, we aim to en- 
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Fig. 2. — Color-magnitude diagram in (i', i'— z') constructed from our deep Suprime-Cam observations. The selection box for follow- 
up spectroscopy is shown with dash-dotted lines. Objects chosen for spectroscopy (plus) and confirmed or possible M-type mem bers 
of NGC 1333 (cros s, see Table [T] and [JJ are marked. Spectroscopically confirmed BDs from the surveys by Wilkins ct al. (2004) and 
IGreissl et all <2007l) are marked with diamond and triangle symbols respectively. 



sure that our candidates show the colors expected for 
reddened very low-mass members of NGC 1333. In Fig. 
[3] we show the color-magnitude diagram constructed by 
combining the Suprime-Cam i'-band photometry with 
the MOIRCS J-band data. The diagram shows all ob- 
jects with data in both bands, but objects selected for 
spectroscopy and the 28 objects confirmed by the spec- 
troscopy (Sect. I33)) are marked. A few of the confirmed 
objects are not in our near-infrared catalogue due to sat- 
uration; for them we use 2MASS data instead. 

For comparison, we overplot the BCAH98, COND03, 
and D USTY00 evol ut ionary t racks for 1 and lOMyr 
from iBaraffe et al.l (|1998t i: IChabrier et all (|2000l ); 
iBaraffe et al.1 ( 2003T ) and the photomet ry for the sub- 
stella r population in the cOri cluster (jCaballero et al.l 
2007). In both cases, we converted the Cousins I-band 
data to the Sloan i '-band using the transformation 
given by Uordi et all ((2006). In principle this requires 
knowledge of the R — I color, which is not always 
available for the crOri sources. We assumed an average 
R — I = 2.0 in these cases; the individual (i'-J) colors 
for a Ori objects may therefore be incorrect by about 
±0.2. Inconsistencies between the Suprime-Cam i'-band 
with the standard bandpass may cause additional 
uncertainties. The a Ori photometry was shifted to the 
distance of NGC 1333, assuming a differential distance 
modulus of 0.7, consistent with recent estimates. aOri 
is a cluster with little extinction; we expect these objects 
to be mostly free from reddening. 



The plot demonstrates that the spectroscopy candi- 
dates cover the color space that is limited on the blue 
side by models and the a Ori sources. This is as ex- 
pected, since the cloud cores in NGC 1333 have signifi- 
cant and variable extinction so the population of young 
cluster members as well as background objects should be 
scattered towards the red side. Apart from reddening, 
there is no systematic offset between our sample and the 
comparison data. Thus, we confirm that our selection 
of candidates covers the color space expected for very 
low-mass members of NGC 1333. 

The intrinsic near-infrared color J — K of late-type 
objects is mostly independent of spectral type and lu- 
minosity; thus it can be used to estimate the extinc- 
tion for the suspected cluster members. This is done 
only for the objects confirmed as having spectral type 
M or later (Sect. I3.3[) . We assume here an intrinsic 
color of J — K = 1.0 mag. Within ±0.2 mag, this is 
consistent with the empirical values for main sequence 
dwar fs and giants with spectral type M (|Bessell fc Brettj 
1988), with the predictions from the evolutionary tracks 
COND03 (for 0.003-0.1 M ), BCAH98 (for 0.03-0.5 Mq), 
and DUSTY00 (for 0.006-0.1 M ), and with the colors of 
brown dwarfs in the 3 Myr old a Ori star- forming region 
down to 0.013 M (jCaballero et al.l I2007D . For young 
objects with masses below the Deuterium burning limit 
there are some indications that the J — K color inc reases 
to 1-1.5 (|Lodieu et alJl2006t ICaballero et al.ll2007f k this 
may lead to an overestimate of the extinction by as much 
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Fig. 3. — Color-magnitude diagram in (J, i'-J) constructed from our MOIRCS and Suprime-Cam observations. The candidates selected 
for spectroscopy from the optical diagram and the confirmed M-type objects (see Table [T] and [2} are marked with plusses and crosses, 
respectively. For com parison, we overplot model isochrones for 1 and 10 Myr and the photometry for the substellar objects in a Ori 
(Caballcro ct al. 2007), shifted to the distance of NGC 1333. For the mo dels and the crOri objects the Cousins-I photometry was shifted 
to the Sloan-i' band using the transformation given bv lJordi et all pOOTJ ). 




J— K (mag) 

Fig. 4. — Near-infrared color-magnitude diagram for the 28 ob- 
jects listed in Table IT1 and l2l The plusses show the MOIRCS pho- 
tometry, crosses are 2MASS values for the cases not covered by our 
data, the circles the dereddened values, assuming intrinsic colors 
of J — K = 1.0 mag (dotted line). The solid lines illustrate the 
reddening path. 

as 2.5 mag for the lowest mass objects. 

For the extinction estimate and throughou t this paper, 
we use the extinction law from lMathisI (|1990f ). with Ry = 
4. In Fig. 2] we plot the (J, J — K) color-magnitude 
diagram with the solid line indicating the reddening path 
and the dotted line showing the assumed intrinsic color. 
For the 28 objects with spectroscopic confirmation (listed 



in Tables Q] and [2]) , we plot our photometry (crosses) or 
2MASS data (plusses), as well as the dereddened J-band 
photometry (circles). For the optical extinction in our 
confirmed sample, we obtain Ay = — 13, i.e. we are 
probing deeper into the cloud than lWilking et ail (2004, 
Ay < 7). The results for Ay are given in Tables [T] and 

m 

This range of optical extinctions is in good agreement 
with the extinction maps as derived from submillimc- 
ter continuum observations (|Enoch et al.l 120061 ). Con- 
verting a 1.1mm map to an extinction map results in 
Ay < 17 mag for 5' resolution for NGC 1333. Our ex- 
tinction estimates for individual objects are all in this 
range, perhaps with the exception of object #12 (see 
discussion in Sect. 13 . 3|) . The comparison also indicates 
that we do not miss a significant number of additional 
objects at even higher extinctions. Note that the maxi- 
mum extinction is lower when derived from 2MASS near- 
infrared data or from 13 CO observations, mainly because 
these indicators do not trace the densest regions of the 
cloud well. 

The main error sources in the extinction estimates are 

a) uncertainties in the estimate of the intrinsic colors and 

b) possible excess emission due to disks and accretion. 
As pointed out above, except for the lowest mass ob- 
jects, the intrinsic colors are likely to be accurate within 
0.2 mag, which causes errors in Ay of < lmag. At the 
lowest masses, the error may b e larger. The disk frac - 
tion in NGC 1333 is high (83%. iGutermuth et al.ll200l . 
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hence a substantial fraction of cluster members may show 
enhanced K-band flux due to excess emission from dust 
in the disk, which leads to a systematic overestimate of 
the extinction. On the other hand, accretion induced 
hot s pots can enhance the J- band as well as the K-band 
flux (jFolha &: Eme rson 199S|). Prior to a more detailed 
characterisation, these effects are difficult to take into ac- 
count. We conservatively estimate that our typical un- 
certainty in Ay is 2 mag; larger errors are possible in 
individual cases. The extinction values will be used in 
the spectral fitting in Sect. 13.31 

3.3. Spectroscopy 

With the multi-object spectroscopy described in Sect. 
12.31 we obtained low-resolution spectra for 53 objects in 
total. The sample of objects with spectra includes 36 
candidates from the primary photometric sample dis- 
cussed in Sect. 13.11 In addition, we have spectra for 
nine objects that are slightly bluer in optical colors than 
our primary candidates sources, to verify the validity of 
the optical color criterion (called 'blue' objects in the fol- 
lowing). Among our spectroscopy sample are nine con- 
firmed young M dw arfs from the sample published by 
iWilkine et aLl(l2?)0l . six of them also in our own photo- 
metric candidate sample. Finally, we included five ran- 
dom objects in our MOS fields. We had 12 additional 
slits on objects that are not well detected and did not 
allow a reliable extraction. All these undetected sources 
have been included as either blue or random objects, and 
do not come from our photometric sample. 

Our goal is to identify young members of NGC 1333 
with substellar masses, i.e. with spectral types later than 
M5. In addition to these young brown dwarfs, our sam- 
ple may contain embedded stellar members of NGC 1333 
and reddened background stars with spectral types ear- 
lier than M5, and (less likely) late M and early L type 
field objects in the foreground. We show some examples 
of the expected shape of the spectra for reddened young 
brown dwarfs, derived from model spectra, in Fig. O 

The low resolution of the spectra does not allow us to 
measure narrow atomic features. Instead, we focus on the 
broadband appearance, i.e. the shape of the H- and K- 
band as well as the relative intensity in these two bands. 
We benefit from the fact that the overall shape of the 
near- infrared spectrum changes distinctly at early/mid 
M spectral types. While earlier type stars show a broad- 
band spectrum which is monotonously decreasing to- 
wards longer wavelength, this is not the case for >M5 
objects, due to the appearance of strong H 2 absorption 
features at 1.3-1.5 /xm and 1.75-2.05 \xm. Most notably, 
this causes the H-band spectrum to have a distinct and 
broad peak centered on 1.68 /im, in contrast to earlier 
type objects, which increases in strength towards L-type 
objects. This fundamental chang e is discussed i n deta il 
in the literature, see for example ICushing et al.l ([2005) . 
As demonstrated in Fig.0 the H-band peak is dominant 
for temperatures lower than 3000 K. While the strength 
and shape of the H-band peak is primarily determined by 
temperature (and gravity, see !4.2j) . there is an additional 
effect of extinction - high extinction (see Sect. 13. 20 . will 
make the H-band peak slightly sharper. As can be seen in 
Fig. [5j this is clearly a minor effect and is only apparent 
when comparing over a wide range of extinctions. 

All spectra were examined visually to look for the typ- 




14000 16000 18000 20000 22000 
Wavelength (A) 

Fi g. 5. — Reddened model spectra from the AMES DUSTY se- 
ries (Allard ct al. 2001) for logg = 3.5, corresponding to pre-main 
sequence ages. Shown are three different temperatures (from top to 
bottom: 3000 K, 2400 K, 1800 K), each with three different extinc- 
tions Ay = 0, 5, 10 mag as solid, dotted, dashed lines. The H-band 
peak and its dependence on temperature is clearly visible. The rel- 
ative strength of H- and K-band is strongly affected by extinction. 
Dash-dotted lines show the spectral range in our observations used 
for classification. 

ical signature of water absorption, mainly the character- 
istic peak in the H-band. The sample clearly falls into 
two groups, where 28/53 objects (53%) show an indi- 
cation of a peak while the others do not. In this clas- 
sification, we erred on the side of caution, and there- 
fore also included objects with only tentative evidence 
of water absorption. This sample of 28 objects is the 
'confirmed' sample of Figs. [20 and is listed in Tables [1] 
and [2] together with the photometric and spectroscopic 
properties. All objects without a clear H-band peak are 
most likely stars with spectral types significantly earlier 
than M5. Our spectra do not permit definitive classi- 
fication of these sources. Out of the 28 confirmed ob- 
jects, 21 come from our photometr i c cand idate sample. 
All nine objects from lWilking et aTl (|2004l ) with spectral 
types M2-M8 are correctly classified as having an H-band 
peak. The 'peak' objects exhibit other spectral features 
of late-type sources, most notably the CO bandheads at 
> 2.29 /im. We do not find evidence for CH4 absorption 
at 1.6-1.8 fim, typical for late L and T dwarfs. 

We compared the 28 spectra with the H-band 'peak' 
with the latest version of th e spectra from the AMES 
DUSTY model atmospheres (|Allard et al.ll2001l ). which 
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Fig. 6. — Spectra of late-type objects in NGC 1333 (solid lines) in comparison with models (dashed lines), first part. The running 
number is the same as in Tables [T] and [2] Shown is the best fit model with the effective temperature as given in the tables. Objects 1-19 
are probable substellar members of NGC 1333. 

are based on the PHOENIX stellar atmosphere code. In 
contrast to the NEXTGEN spectra from the same code, 
DUSTY spectra include dust formation in the convective 
envelope and dust opacities. The more recent DUSTY 
models also include more complete H2O and TiO line 
lists. To probe if the specific treatment of dust is rele- 



vant, we compare d the DUSTY spectra with the DRIFT 
PHOENIX series (|Schmidt et al.ll2008f) . which include a 
dust model and consistently treat cloud formation and 
its feedback on the atmosphere. For the resolution and 
wavelength coverage comparable with our observed spec- 
tra, there is no significant difference in the temperature 
regime > 2200 K. For lower temperatures, the dust treat- 
ment beco mes critical and alter native models need to be 
considered (|Helling et al.ll2008l ). 

We used the DUSTY model spectra with temperatures 
of 1800-3900 K in steps of 100 K and the lowest available 
log<? (3.5 or 4.0). The model spectra were binned and 
smoothed in the same way as the science spectra. All 
spectra are scaled to a flux of 1.0 at 1.7 fim. Each ob- 
ject spectrum was then compared with the model series, 
reddened with the optical extinction derived in Sect. 13.21 
The quality of the match was judged by calculating the 
X 2 aided by visual inspection. To avoid being misled by 
excess continuum due to disks or accretion and uncer- 
tainties in the extinction, we put emphasis on matching 
the H-band shape, rather than the overall spectrum. 



In Figs. [6] and [7] we show all 28 object spectra with 
H-band peak (solid lines) and their best fit model spec- 
trum (dashed lines). The ID numbers are the same as in 
Table [1] and O The best match temperatures are given 
in the tables. For 18 out of 28 objects, this procedure 
yields a minimum x 2 of < 0.5 and a well-matched spec- 
trum, without further adjustments. Four more objects 
(#4, #10, #20, #24) fulfill the same criteria after minor 
changes to the extinction value (A Ay < 2 mag). Object 
#26 requires increasing the extinction from Ay = 4.4 
to 8 mag, to obtain a good match, which results in 
T = 3800 K. As argued in Sect.OJ the status of this ob- 
ject is unclear; it may be a red giant in the background 
of the cluster. Object #12 shows steeply increasing flux 
levels in the K-band and all signs of strong reddening. 
We obtain a good match between model and observed 
spectrum by increasing the extinction to Ay = 18 mag. 
The source is only partially covered in the H-band, there- 
fore the temperature estimate is preliminary. Based on 
the slope of the H-band spectrum, the best estimate is 
T = 2200 - 2800 K. Two objects (#8, #18) agree well 
with the model in the H-band, but show continuous ex- 
cess emission in K-band with respect to the best fitting 
model. Both sources are found to have mid-infrared color 
excess due to the presence of circumstellar material (Sect. 
14. 5| . which could explain the K-band excess. 

Thus, for 26 out of 28 objects we obtain a decent match 
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Fig. 7. — Spectra of late-type objects in NGC 1333 (solid lines) in comparison with models (dashed lines), second part. The running 
number is the same as in Tables [T] and [2] Shown is the best fit model with the effective temperature as given in the tables. Objects 1-19 
are probable substellar members of NGC 1333. 



between the model and observed spectra. In Fig. [8] the 
typical uncertainty of the fitting procedure is illustrated 
by comparing the observed spectra for three objects from 
Table [T] with models for various effective temperatures. 
The x 2 values used to select the best-fitting model are 
listed as well. As can be seen in the figure, the procedure 
often does not allow us to distinguish reliably between 
two models with AT = 100 K. However, for AT > 200 K 
there are discernible differences, particularly in the ri- 
band peak, between model and observation, combined 
with a clear increase in \ 2 ■ Thus, a conservative estimate 
for the typical uncertainty in T e g is ±200 K. For all >M5 
object s, the spectral types - either from Wilking et al. 
(2004) or from the spectral index, see below - are con- 
sistent with the temperatures from model fitting within 
±200 K or ±1 s ubtyp e (using the temperature scale from 
iLuhman et all (|2003h ). Temperatures of > 3500 K have 
to be treated with caution, as our model grid is limited 
at 3900 K. 

The two objects with unsatisfactory results in the fit 
arc #11 and #15. For #15, the H-band spectrum is 
well- matched with a temperature of 3000 K which is con- 
sistent with the spectral type: M6.5 from Wilking et al. 
(|2004ft and M7.5 from our spectrum, see below. The 
literature K-band spectrum shows a flux increase from 
2.1 to 2.3 /i 77i, as expected for a mid/late M-type object. 
Thus, we consider the temperature estimate of ~ 3000 K 



to be reliable. Our K-band spectrum declines with wave- 
length, indicative of a significantly higher temperature, 
which prevents us from obtaining a good fit. The object 
also has a mid-infrared color excess, indicative of circum- 
stellar material (Sect. |4~5|) . Source #11 exhibits strong 
emission at 2.12 /im, consistent with the 1-0 S(l) line of 
molecular hydrogen at 2.122 /ito, seen in protostellar out- 
flows. At A > 2.25 fj,m, the flux decreases sharply, which 
is not reproduced by any models or literature spectra in 
this temperature domain. Based on the H-band peak, 
our best estimate for the effective temperature is 2600 K, 
but this has to be treated with caution. 

As a complementary approach, we c ompu ted the spec- 
tral index proposed by lAllers et al.1 (|2007f ) for spectral 
classification of young BDs, applicable from M5 to L5. 
This index has two distinct advantages compared with 
the variety of other index schemes in the literature, as 
it is gravity insensitive and does not require wavelength 
coverage in regions heavily affected by telluric water ab- 
sorption. The spectr al types from this index confirm the 
Wilki ng" et al.1 (|2004h spectral types for the seven >M5 
objects within ±1.1 subtypes. Where applicable, we use 
this index as an additional indicator for the spectral type 
of our objects (see Table [1} . 

3.4. Summary of substellar member selection 
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Fig. 8. — Illustration of the spectral fitting process: Each panel shows the spectrum of a probable brown dwarf in NGC 1333 from 
Table [T] in comparison with model spectra for a range of effective temperatures. For clarity, the models are offset by 0.05 units. For each 
comparison, we give the effective temperature of the model and the \ of the fit. From the \ 2 values we infer a conservative uncertainty 
of ±200 K for the effective temperature. The plot demonstrates that larger offsets in T e ff lead to visible discrepancies between model and 
observation, particularly in the shape of the H-band peak. Note: For object #4, the photometrically determined extinction of Ay = 2.3 mag 
gives x 2 •> 0.37; much better results are achieved with a slightly lower Ay. We use here Ay = 1 mag. 



In summary, from our 36 photometrically selected can- 
didates, we confirm 21 as M-type objects. In addition, 
three out of nine bluer objects and one out of five ran- 
dom objects are cl assified as M-type sou rces. Nine clus- 
ter members from IWilking et al.l (|2004f ) are confirmed 
by our spectroscopy, seven of them overlapping with our 
photometrically selected sample. In total, we have 28 
confirmed objects with a H-band peak, indicative of spec- 
tral type M. In this confirmed sample are 19 objects with 
spectral types >M5 or effective temperatures of <3000K, 
which would qualify them as brown dwarfs, if they are 
in fact young members of NGC 1333 (see Sect. 14. 2p . 
This subsample is summarized in Tabled the remaining 
nine objects in Tabled For completeness , we note that 
four p hotometric candidates reported by IWilking et al.l 
(2004) do not show the H-band peak, i.e. are not clas- 
sified as M-type sources: MBO 39, 59, 76, 109. These 
objects may still be cluster members with earlier spectral 
types. 

4. CHARACTERIZATION OF CONFIRMED LATE-TYPE 
OBJECTS 

4.1. Effective temperatures 

We have determined effective temperatures for the 
sample of late-type objects in Tables Q] and [2] by fit- 
ting the observed to model spectra. Out of 28 sources, 
19 have temperatures between 2500 and 3000 K, for the 
other nine we obtain higher values. In addition, we have 
already derived extinctions, based on the J — K colors, 
and dereddened J-band magnitudes (see Sect. 13.21 Fig. 

a. 

We use the dereddened J-band photometry to obtain 
an alternative estimate of effective temperatures by com- 
paring with the BC AH98 evolutionary tracks BCAH98 
(|Baraffe et al.1 11998) for the typical cluster age of 1 Myr. 
This relies on the assumption that the objects are in 



fact young members of the cluster NGC 1333, for which 
we adopt a distance of 300 pc, c o nsistent with the es - 
timates bv lde Zeeuw et al.l (| 19991 ): iBelikov et all (|2002h . 
The J-band is the wavelength regime that is normally 
assumed to be least affected by accretion, disk excess, 
and magnetic activity. It is found that the BCAH98 
temperature vs. J-mag relation is well approximated by 
a fourth degree polynomial. Applying this function to 
the dereddened J-band photometry yields 'photometric' 
temperatures from 2100 to 3400 K. 

For 15 objects, the photometric and spectroscopic tem- 
peratures agree within ±200 K, confirming our assess- 
ment of T e ff given in Sect. 13.31 Four objects (23, 26-28 
in Table [2|) have photometric temperatures 800-1200 K 
lower than the ones from the spectroscopy, which can be 
explained with being background objects. Assuming in- 
trinsic colors for dwarf stars from the 1 Gyr isochrone in 
the BCAH98 model tracks puts their distances at 360- 
950 pc. The remaining objects have offsets of 200-500 K 
between the two values, and in most cases the photomet- 
ric temperatures are too low. One possibl e explanation is 
an age spread in the cluster. As shown by I Winston et all 
(|2009l) . the members of NGC 1333 are significantly dis- 
persed in the HR-diagram, possible indicating an age 
spread from < lMyr to several Myrs. Furthermore, the 
model tracks have wel l-documented uncer tainties in this 
mass and age regime (TBaraffe et al.ll2002Th 

For two objects, #23 and #26, there is a clear discrep- 
ancy between our spectroscopic temperatures of 3800 and 
39 00 K and the spectral types of M2.6 and M4.2 derived 
by IWilking et alJ (|2004l ). which would imply tempera- 
tures of only 3500 and 3300 K, respectively. As pointed 
out above, these are likely to be backgroun d objects. The 
K-ban d index for spectral typing used by Wilking ct al. 
(2004) gives significantly earlier spectral types for giants 
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TABLE 1 

Probable substellar members of NGC 1333. 



HO. 




(J2000) 


<5(J2000) 


i ' ( mfif] 


7, 3 ( Tna en 


.T (iTlflP') 




Av 


+ crt 


SpT b 


Notes c 


1 


03 


28 


47.66 


+31 


21 


54.6 


23.540 


21.343 


17.55 


15.24 


6.9 


2800 




MB0139 


2 


03 


28 


54.93 


+31 


15 


29.1 


21.040 


18.922 


15.99 


14.22 


4.1 


2600 


M6.5 


cx d > c 


3 


03 


28 


55.25 


+31 


17 


35.4 


19.386 


17.602 


15.09 


13.43 


3.5 


2900 




ASR38 d 


1 


03 


28 


56.50 


+31 


16 


03.1 


22.999 


21.233 


18.17 


16.73 


2.3 


2500 




e 


5 


03 


28 


56.94 


+31 


20 


48.6 


19.417 


17.935 


15.48 


13.92 


3.0 


2900 


M6 


MB091, ex 


6 


03 


28 


57.11 


+31 


19 


12.0 


21.157 


19.545 


17.24 


15.34 


4.8 


2700 


M8 


MB0148, ex 


7 


03 


28 


58.43 


+31 


22 


56.7 


19.230 


17.856 


15.21 


13.54 


3.5 


2800 


M6.5 


MBO80, ex 


8 


03 


29 


03.39 


+31 


18 


39.9 


20.059 


18.199 


15.85 


14.03 


4.3 


2600 


M8.5 


MBO88, ex 


9 


03 


29 


05.56 


+31 


10 


13.7 


21.288 


19.425 


16.91 


15.76 


0.8 


2600 


M8 




10 


03 


29 


05.64 


+31 


20 


10.7 


20.080 


18.842 


17.11 


15.49 


3.3 


2500 




MB0143, ex d ' c 


11 


03 


29 


07.17 


+31 


23 


22.9 


22.989 


21.187 


17.87 


15.62 


6.6 


(2600) 




MB0141 e 


12 


03 


29 


09.33 


+31 


21 


04.1 


22.553 


20.336 


16.33 


13.29 


10.7 


(2500) 




MBO70, ex c 


13 


03 


29 


10.79 


+31 


22 


30.1 


19.483 


17.778 


14.96 


13.05 


4.8 


3000 


M7.5 


MB062 


14 


03 


29 


14.43 


+31 


22 


36.2 


18.104 


16.960 


14.55 


13.02 


2.8 


2900 


M7 


MBO66 


15 


03 


29 


17.75 


+31 


19 


48.1 


19.091 


17.326 


14.80 


12.99 


4.3 


3000 


M7.5 


MB064, ex d ' c 


16 


03 


29 


28.16 


+31 


16 


28.4 


16.934 


15.808 


13.22 


12.53 


0.0 


2600 


M7.5 


random 


17 


03 


29 


33.88 


+31 


20 


36.1 


20.383 


18.848 


16.47 


15.37 


0.5 


2500 


M8 


MBO140 


18 


03 


29 


35.71 


+31 


21 


08.5 


22.633 


21.089 


18.50 


16.94 


3.0 


2500 




blue, ex 


19 


03 


29 


36.36 


+31 


17 


49.8 


22.528 


20.716 


17.91 


16.38 


2.8 


2700 







a Estimated by comparing the spectra to models. Uncertain fit results are in brackets. 

Estimated using the spectral index and calibration from lAllers et ~ ] d2007lh if appli cable. 

c 'ex' means mid-infrared color excess, ASR and MBO object id's arc candidates from Wilking ct al. (2004). 
^ Near-infrared photometry from 2MASS. 

c Temperature estimate changed after visual inspection, sec Sect. HT3l 



TABLE 2 

Possible stellar members of NGC 1333. 



no. 


a(J2000) 


<5(J2000) i 


i' (mag) 


z' (mag) 


J (mag) 


K (mag) 


Ay 


Trff a SpT b 


Notes c 


20 


03 28 43.23 


+31 10 42.6 


21.864 


19.606 


15.37 


12.03 


12.3 


3400 


ex d > e 


21 


03 28 47.34 


+31 11 29.9 


21.626 


19.330 


15.48 


12.70 


9.4 


3100 


d 


22 


03 28 50.97 


+31 23 47.9 


20.380 


18.418 


14.66 


11.24 


12.7 


3500 


MB026 d 


23 


03 29 02.16 


+31 16 11.1 


17.197 


16.245 


14.53 


13.57 


0.0 


3900 


ASR03, giant? d 


24 


03 29 08.32 


+31 20 20.3 


19.376 


17.809 


14.80 


12.48 


7.0 


3100 


MB044 C 


25 


03 29 11.64 


+31 20 37.5 


20.544 


18.738 


15.43 


12.67 


9.3 


3400 


MB058 d 


26 


03 29 16.75 


+31 23 25.3 






15.42 


13.58 


4.4 


3800 


MB079, ex, giant? d ' c 


27 


03 29 22.36 


+31 21 36.8 


21.407 


20.197 


17.71 


15.89 


4.3 


3300 


MB0177, blue 


28 


03 29 33.50 


+31 17 56.5 


22.926 


21.477 


18.52 


16.27 


6.6 


3700 


blue 



a Estimated by comparing the spectra to models. Uncertain fit results are in brackets. 

b Estimated using the spectral index and calibration from lAllers et all i2007f) , if applicable. 

c 'ex' means mid-infrared color excess, ASR and MBO object id's are candidates from Wilking ct al. (2004). 

d Near-infrared photometry from 2MASS. 

c Temperature estimate changed after visual inspection, sec Sect. [373] 

than for dwarfs (about M0 or earlier for these two ob- 
jects). These types would be consistent with our effec- 
tive temperatures. Thus, these two objects could be red 
giants in the background of NGC 1333. However, assum- 
ing the intrinsic giant colors from Bessell fc Brettl £1988), 
this would put them as an unreasonably large distance 
of 30kpc for spectral type MO. A significantly earlier 
spectral type is at odds with our spectroscopy. Thus, 
the nature of these objects remains uncertain. 

4.2. Cluster membership 

For the likely brown dwarfs in our sample (see Table 
[T]), the best evidence for youth and thus cluster mem- 
bership is given by the shape of the H-band peak in 
the spectra. For substellar objects, the H-band peak 
appears to be gravity sensitive. Low gravity and thus 
youth is signified by a distinctly sharp peak, as op- 
posed to the more rounded peaks seen in old objects 



(jKirkpatrick et alJ l2006). The same characteristic is also 
seen in the available model spectra (jAllard et al.ll200"il ). 

All 19 sources with effective temperatures of 3000 K 
or lower listed in Table [T] show a sharp peak. As an 
example, we show four of them in comparison with pub- 
lished spectra for field dwarfs of similar spectral type 
from the NIRS PEC Brown Dwarf Spectroscopic Survey 
(|McLean et al.ll2003l) in Fig. [SJ left panel. The difference 
in the shape of the H-band feature is obvious. Particu- 
larly, the slope on the red side of the peak is significantly 
steeper in our objects. On the other hand, our spectra 
agree well with those of young b rown dwarfs in Taurus, 
taken from iMuench et all (|2007t ). as shown in Fig. [9l 
right panel. 

Statistically, we do not expect a significant contami- 
nation with late field M dwarfs in our sample in Table 
[1] Ac cording to the data compiled in ICaballero et al.l 
(2008), the space density of M6-7 dwarfs is 4.9-10" 3 pc _1 
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and drops steeply to 1.3- 10 3 pc 1 for M8-9. Our sample 
covers an area of about 200arcmin 2 which corresponds 
to ~ 1.5 pc 2 at the distance of NGC 1333. In a cone with 
this footprint and 300 pc height we expect about one M6- 
M9 field dwarf. Thus, for the M6 or later objects in Table 
Q]the contamination by field objects is negligible. 

We conclude that the >M5 objects identified in this 
survey ('probable' members in Table [JJ are in fact young 
substellar members of the cluster NGC 1333. For the 
earlier type stars ('possible members' in Table 0, sig- 
nificant contamination by field dwarfs or giants is more 
likely, as their space density is larger and the H-band 
signature of youth is less pronounced. Further evidence 
for youth comes from the fact that many of our objects 
show mid-infrared color excess indicative of the presence 
of circumstellar material, as discussed in Sect. 14.51 and 
from the spatial clustering around the main cloud core 
in NGC 1333, see Sect, 

4.3. Masses 

We estimate the masses of the probable brown dwarfs 
in our sample by comparing the effective temperatures 
and the dereddened J-band photometry with the theo- 
retical model tracks for an assumed age. Considering the 
uncertainties in model tracks and object ages, we refrain 
from deriving masses for individual objects. Instead, we 
focus on establishing the approximate mass limits in our 
sample. 

The range of effective temperatures in this sample is 
2500-3000 K. Comparing directly to the 1-3 Myr tracks 
from BCAH98 yields masses of 0.02-0.1 M . The result 
is very similar for the COND03 and DUSTY00 tracks. If 
some objects are younger than 1 Myr, they will be slightly 
more massive for a given temperature. Taking into ac- 
count 200 K uncertainty in our temperature estimate, the 
lowest mass limit becomes 0.012-0.015 Mq, according to 
the COND03 and DUSTY00 tracks. Based on the dered- 
dened J-band photometry, the mass range in our brown 
dwarf sample is 0.006-0.1 Mq for the 1 Myr isochrones, or 
0.012-0.2 Mq for the 5 Myr isochrones. Note that the lat- 
ter values are not significantly affected by uncertainties in 
the distance to NGC 1333. The most recent distance esti- 
mates agree f airly well at ~ 300 pc with an uncertainty o f 
about 30pc (|de Zeeuw et al .1119991 : iBelikov et all 12002( 1 . 
corresponding to ±0.2 mag or about 10% in object mass. 

There are two reasons why we put more trust in the 
mass estimates obtained from the temperatures: 1) Lu- 
minosities depend much more on age than temperatures 
on the pre-main sequence, thus the mass estimates based 
on temperatures are more robust again st age spread, 
which may be significant in this cluster ((Winston et alJ 
2009). The Lyon evolutionary tracks are known to un- 
derestimate object masses derived from fluxes. For ex- 
ample, foj_th^_j^ojong_brown dwarf eclipsing binary in 
ONC ((Stassun et alJ 12006( 1 the BCAH98 tracks predict 
a mass of 0.055-0.075 Mq based on its J-band magni- 
tude (assuming a distance of 435 ± 55 pc); the actual 
system mass is 0.09 Mq. More indirectly. IMohant y et al.l 
(2004) find that the models predict brown dwarfs to be 
too bright in J-band at a given mass, by up to one or- 
der of magnitude. Both tests indicate that the masses 
from the model tracks, when inferred from luminosities, 
may systematically underestimate the actual masses by 
a factor of up to two. 



Based on the comparison of effective temperatures with 
the three model tracks, we conclude that the lowest mass 
limit in our brown dwarf sample is 0.012-0.02 Mq. At the 
high mass end, the mass range of the objects in Table [1] 
extends across the substellar boundary into the stellar 
regime. Thus, a few of these objects may turn out to be 
very low-mass stars, instead of brown dwarfs. 

4.4. Spatial distribution and binarity 

In Fig. [10] we show the position of our photometric can- 
didates and the confirmed members of NGC 1333 in rela- 
tion to the gas in the cloud and other samples of cluster 
members. As already pointed out in Sect. 13. 1( the photo- 
metric candidates show a strong clustering towards the 
center of the cloud. On one side, this can be explained 
with the increasing extinction, causing background ob- 
jects to appear heavily reddened. On the other hand, as 
seen in the plot, YS Os candidates selected fr om Spitzer 
mid-infrared excess (|Gutermuth et alJ l2008f) show the 
same strong clustering, indicating that the distribution 
of members is very compact and does not extend to the 
edges of our survey field. Our spatial coverage thus fully 
encompasses the extent of the cluster. The small num- 
ber of photometric candidates in the outer region is also 
evidence for the lack of contamination in our sample, as 
argued in Sect. 14.21 

All our spectroscopy fields were located within 
±0.1 arcmin of the cluster center, an area which contains 
the majority of known cluster members. As discussed in 
Sect. 13.31 the 'success rate' in finding late-type sources is 
58% (21/36) in the photometrically selected sample, in- 
dicating a substantial rate of contamination among the 
remaining photometric candidates. Due to our limited 
spatial coverage in the spectroscopy, we cannot give con- 
straints on the nature of the (few) candidates outside the 
cloud core. Since the surface density of members drops 
rapidly with increasing distance from the center, we ex- 
pect most of the outlying candidates to be foreground or 
background sources. 

The very compact distribution of brown dwarfs agrees 
with the more detailed analys is for the distribution o f 
stars in NGC 1333 provided bv lGutermuth et all ((2008(1 . 
who report a steep decline of the density of members for 
cluster radii larg e r than 0.3 pc (0.06 degree at 300 pc). 
iGutermuth et all (|2008( l also identify clear asymmetry 
in the distribution, attributed to the structure of the 
underlying cloud. This all suggests that the cluster has 
experienced only limited dynamical evolution and most 
of the population in the core is in fact very young. Ages 
ranging from less than 1 Myr to 3 Myr have been inferred 
from HR diagram fitting for the objects conc entrated in 
the molecular cloud gas ((Winston et alJ[2009f ). 

We visually checked the deep i'- and z'-band images 
for companions within a radius of 3.0" of each of our 
spectroscopically confirmed M-type sources. The lower 
separation limit for the detection of companions is set 
by the seeing and the rejection criterion for elongated 
objects (see Sect. 12. ip : we determine this to be at 
~1.0" (Sect. [23]) . corresponding to a separation range of 
300 to 1000 AU for the distance of NGC 1333. No com- 
panions were found. For substellar objects in NGC 1333 
(Table [TJ this gives a 2a upper limit for binarity in the 
given separation range of 19%. 

A number of recent studies clearly show a dearth of 
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Fig. 9. — Spectra of very low-mass objects in NGC 1333 from Tablep] (solid) in comparison with spectra of field dwarfs from the NIRSPEC 
Brown Dwarf Spectroscopic Survey (left panel, McLean ct al. 2003) and with young brown dwarfs in Taurus (right panel, IMuench et al.l 
l2007t) . The field dwarfs are from bottom to top Wolf 359 (M6), VB 10 (M8), LHS2065 (M9); the Taurus BDs are from bottom to top 
MHO 4 (M7), KPNO 1 (M7.5), and KPNO 5 (M8.5). The objects in NGC 1333 have a significantly 'sharper' peak than the field dwarfs, but 
are well matched by the young Taurus BDs. To account for the higher extinctions in NGC 1333, the literature spectra have been reddened 
to Ay = 5 mag. 



wide binaries in the very low-mass regime. In the central 
parts of the Cha I star forming region, for example, a 
binary fraction of 1 1% is found for separations > 20 AU 
(Ahm ic~et al.| [2007), and zero systems with separations 
larger than 50 AU. Our result in NGC 1333 is consistent 
with these findings ( for a review on VLM binarity see 
iBurgasser et alJl2007f >. 

4.5. Spitzer photometry 

The NGC 1333 cluster is part of the Perseus star- 
forming complex and thus has been covered by the 
Spitzer Legacy Program 'From Cores to Disks' (C2D, PI: 
Neil Evans). From the C2D source catalogues (version 
S13, full Perseus catalogue) we obtained the IRAC and 
MIPS fluxes from 3.6 to 24 [im for all objects in Tables 
[T]and [2] Out of 28 sources, 16 are well-detected with er- 
rorbars < 0.2 mag in all four IRAC bands. The standard 
IRAC color-color diagram is shown in Fig. [TTJ The color 
range for Class II objects is indicated. Eleven objects 
clearly show IRAC color excess, and thus are likely to be 
surrounded by circumstcllar material. Ten of them are 
also safely detected at 24 /im, with magnitudes > 2 mag 
larger than at 8 \xm. This confirms youth and thus mem- 
bership to NGC 1333 for a large fraction of our sample. 
Objects with mid-infrared color excess are marked with 



'ex' in Table [TJ The disk fraction in our (small) sam- 
ple is ~ 70%, consistent with previous findings for stars 
(83%. IGutermuth et all 120081 ) and brown dwarfs (74%, 
IWilking et aHl200l . 

The C2D data h as rec ently been analysed in detail by 
IGutermuth et al.l {2008). In this survey, all 11 objects 
with disks in our sample are classified as Class II sources. 
Two more objects in our sample (# 23 and #25) are 
Class II objects with disks according to IGutermuth et ail 
(2008), but their IRAC errors are too large to be consid- 
ered for Fig. [TTJ 

5. THE MASS FUNCTION IN NGC 1333 

5.1. Ratio of stellar to substellar members 

Com bini ng our survey with t he ones by IWilking et al.l 
(|2004h and iGreissl et all (|2007f) there is now a substan- 
tial sample of spectroscopically confirmed brown dwarfs 
in NGC 1333. We use this to improve the constraints 
on the ratio of low-mass stars to brown dwarfs, a num- 
ber recently measured for a range of other star-forming 
regions. We identify 19 objects as young M6 or later 
so urces. Seven of them have previously been confirmed 
bv IWilking et alj (|2004l ). Their survey has found nine 
more objects with M6 or later spectral type, which we 
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FlG. 10. — Spatial cove rage of our survey and distribution of objects. Crosses are all candidates selected from optical phot ometry with 
i' > 18.5 mag (Sect. 13.11 1. diamonds the spectroscopically confirmed objects, triangles the confirmed bro wn dwarfs fro m Wilkin g et al.1 
(2004). The plusses are YSO candidates selected from Spitzer data without spectroscopic confirmation from iGutcrmuth et al. (2008]). The 
underlying contours show the distribution of 13 CO from the J = 1 — map by Ridge et al. (200rJ). The dotted and dashed lines show the 
area covered in our optical and near-infrared survey, respectively. 




0.5 

[6.5] - [8.0] (mag) 

FlG. 11. — IRAC color-color diagram for the 16 objects from 
Tables [l] and [2] which have errorbars < 0.2 mag in all four bands. 
The approximate locus of the Cla ss II objects is ovcrplotted as 
a dashed square IjAllen et al.l 120041 ). Typical photospheric colors 
for late-type objects are around zero. Eleven objects have colors 
consistent with color excess due to circumstellar material. 

consider to be good candida tes for brown dwarfs. The 
study bv iGreissl et al.l ()2007l ) has limited spatial cover- 
age (3arcmin 2 ), but finds five more objects in this spec- 
tral regime. All these objects are located close to the 
cluster center, which points to a high likelihood of being 
members. Thus, the current total census of likely brown 
dwarfs in this cluster is 33. As there are a number of 
additional photometric candidates in the cluster region 



(see Fig. [10|) that lack spectroscopic confirmation, this 
number should be treated as a lower limit; on the other 
hand, as pointed out in Sect. 14. 3[ some of these objects 
may be slightly above the substellar boundary. 

To constrain the number of stellar members in 
NGC 1333, we use t he Spitzer selected sample by 
iGutermuth et alj (|2008fl . which identifies 137 likely YSOs 
in the cluster, 93 of them also have 2MASS near-infrared 
magnitudes. From these, 91 are in our survey area. We 
dereddened the photometry for these objects based on 
the J — K color, using the procedure discussed in Sect. 
13.21 Comparing the dereddened J-band magnitude with 
the BCAH98 evolutionary track allows us to select ob- 
jects in a given mass range. In this sample, 38 sources 
are likely to have masses between 0.08 and 1 Base d 
on the spectroscopic follow-up by iWinston et al.l (2009) , 
this sample is unlikely to have significant field star con- 
tamination. Spitzer only finds objects with disks; in fact, 
the cluster is known to h arbour a significant p opulation 
of disk-less T Tauri stars (IWinston et al.ll2009ft. A ssum- 
ing a disk fraction of 83% ( Gutermuth et al.ll2008f ) gives 
a total number of ~ 50 stellar cluster members in this 
mass regime. Taking into account the uncertainties in 
our mass cuts, this number may be inaccurate by up to 
20%. 

Although the Spitzer sample of stars has been identi- 
fied in a way not comparable to our own brown dwarf 
survey (and the two other ones in the literature), we are 
confident that it gives a robust estimate for the number 
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of stars in this mass range in the cluster. The range of 
extinctions and the spatial distribution are similar in the 
Spitzer sample and the brown dwarf sample. Moreover, 
the disk fraction among the brown dwarfs is similar to 
the one seen in stars (Sect. I4.5[) . Thus, in terms of evo- 
lutionary state and spatial coverage the stellar sample is 
representative. 

From these numbers, we estimate the ratio of (low- 
mass) stars (0.08 < M < 1.0 Mq) to brown dwarfs 
(0.015 < M < 0.08 M Q ) in NGC 1333 to be 50/33, 
i.e. 1.5 ± 0.3. This can directly be compared with 
the literature. In an IMF analysis for seven differ- 
ent star-formin g regi ons (not including NGC 1333), 
lAndersen et all (|2008D find this ratio to be 3.3-8.5 de- 
pending on the region, with the lowest values measured in 
the ONC, NGC 2024, and Chamaeleon. Taking into ac- 
count the la uncertainty in their values, the minimum ra- 
tio becomes 1.9 in Chamaeleon. These numbers are still 
significantly higher than the ratio in NGC 1333, pointing 
to an overabundance of brown dwarfs in NGC 1333 by a 
factor of 2-5. 

This result i s consistent wi t h the previous findings for 
this cluster. iWilking et al.1 (|2004l ) derived a ratio of 
stars (0.08 < M < 1.0 Mq) to very low-mass objects 
(0.04 < M < 0.1 M Q ) of 1.1±8;|, and considered this to 
be a lower limit, i.e. fully comparable with our result. 
(This would become 2.0, if they assume an age of 1 Myr 
instead of 0.5 Myr.) Includ ing their lower mass brown 
dwarfs, iGreissl et all (|2007l ) find a value of 0.64 ± 0.38 
for this ratio. Although these numbers are not fully com- 
parable with our estimate, due to the slightly different 
mass ranges, they are all clearly lower than the respec- 
tive ratios in other regions. Thus, the mass function of 
NGC 1333 has a clear anomaly with an overabundance 
of brown dwarfs relative to low-mass stars. 

5.2. The minimum mass in NGC 1333 

The depth and large spatial coverage of our survey 
allows us to put meaningful limits on the abundance 
of planetary-mass objects in NGC 1333 in comparison 
with other regions. Recapitulating, the coolest objects 
in NGC 1333 have effective temperatures of 2500 K, cor- 
responding to masses of 0.012-0.02 Mq. That means that 
there is no planetary-mass object in our spectroscopically 
confirmed sample. In the following, we will compare with 
the mass completeness limits in our survey and the fre- 
quency of planemos in other star-forming regions. Given 
the considerable uncertainties in the model evolutionary 
tracks in this age and mass regime, we will follow two 
guidelines: a) Conversion from effective temperatures to 
masses is more reliable than from luminosities (or mag- 
nitudes) to masses, see the discussion in Sect. 14.31 b) 
We prefer to base the arguments on relative comparisons 
consistently using the same model tracks, and do not 
place strong emphasis on the absolute mass values. 

Our coolest and lowest mass objects in the spectro- 
scopic sample have T c ff = 2500 ± 200 K, i.e. we rule 
out the presence of objects with 2200 K or cooler in 
this group of objects. The five faintest confirmed brown 
dwarfs have 22.5-23.5 mag in the i'-band with extinctions 
Av ^ 10 mag. Our completeness limit is 24.7 mag in 
i'-band, i.e. 1.2-2.2 mag deeper than the faintest con- 
firmed objects. Based on the lMyr isochrones from 
COND03 and DUSTY00, an object that much fainter 



than one with 2500 K has 0.006-0.01 Mq, with T eff of 
2000-2200 K, i.e. for A v < 10 mag we are complete 
down to ~ 0.008 Mq. For A v < 5 mag the limit be- 
comes ~ 0.004 Af or 1800 K. Thus, our survey is com- 
plete to masses significantly below the low-mass limit 
of the confirmed objects. In fact, six objects with i'- 
mag > 23.5 mag have been observed spectroscopically 
and turned out not to have an H-band peak, indicative 
of low temperatures. 

As discussed in Sect. [TJ previous surveys with compa- 
rable depth (in terms of masses) have been carried out in 
the dOri cluster, the Orion Nebula Cluster (ONC), and 
the Chamaeleon I star-forming region (Cha I). The most 
recent ce nsus of the subste l lar po pulation in a Ori, pub- 
lished by ICaballero et all (|2007f) . lists 34 objects with 
0.015 < M < 0.08 M Q and 12 with lower masses. In 
their mass scale, the planetary-mass objects have spec- 
tral types M9 or later, indicating that they exhibit lower 
temperatures (and thus lower masses) than the coolest 
objects in NGC 1333. The numbers of BDs and planemos 
in a Ori are robust against the particular choice of evo- 
lutionary tracks and un certainties in ag e , dist ance, and 
reddening. According to lCaballero et all (|2007l ). the cen- 
sus is expected to contain ~ 6 contaminating foreground 
dwarfs, ~ 4 of them in the planemo regime. Thus, the 
best estimate for the planemo vs. BD ratio in a Ori is 
8/32 or 25±g%. This is a lower limit, as the survey 
does not show any cut-off of the mass function down 
to 0.006 Mq, their completeness limit. According to re- 
cent adaptive opt ics observatio ns in the core o f a Or i 
(|Bouv et alJ 120091 ) the survey by ICaballero et all (|2007l ) 
might miss a number of objects with planetary masses 
close to the bright central star. 

In the ONC, the highly variable extinction makes an 
unbiased census difficult. In the largest spectroscopic 
survey in the substellar regime to date, [Weights et al.l 
(2009) have confirmed ~ 38 objects with masses at or be- 
low the Hydrogen burning limit. According to their mass 
estimates, 17 objects in this sample have 0.015 < M < 
0.08 Mq and 20 objects have even lower masses. The 
masses are based on a comparison of the observed HR di- 
agram with evolutionary tracks (including COND03 and 
BCAH98, as used in this paper) and thus rely partly on 
spectroscopic temperatures, comparable to our estimates 
in NGC 1333. Deriving the mass limits fr om the dered- 
dened H-band photometry, as given by Weights et al. 
(f2009h . yields 15 objects with 0.015 < M < 0.08 M and 
14 below (assuming a distance of 450 pc). Their analy- 
sis of the mass function thus supports the presence of a 
planetary-mass population in the ONC. 

There are two alternative estimates for the frequency 
of planemos in the ONC in the literature: According to 
iLucas et al.l (|2006f ). 7.5% (1-14%) of the total population 
in the cluster have a mass of 0.003-0.015 M , where the 
wide range of p ossible values is m ostly a reflection for 
age uncertainty. ILucas et al.l (|2005f ) find a drop by factor 
two in the mass function at the Deuterium burning limit, 
a deviation from a flat mass function in the substellar 
regime. Their upper limit for the fraction of planemos in 
the total cluster population is 10-13%. 

For Cha I, t he most up-to-date unbiased census by 
iLuhmanl (2007) yields a total cluster population of 226 
objects, among them 28 brown dwarfs, and four sources 
with masses below the Deuterium burning limit. Based 
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on these numbers, the planemo fraction is 2% relative to 
the total population and the planemo vs. brown dwarf 
ratio is 17%. In this region, however, the census is com- 
plete only down to ~ 0.01 M©. The mass function is 
flat in the substellar regime and does not sho w indica- 
tions for having reached the minimum mass (jLuhmarJ 
2007). Thus, the total number of planemos may be twice 
as high, depending on the cut-off of the IMF. In fact, 
iLuhman fc Muenchl ([20081 ) published one more source 
in the planetary-mass regime based on a Spitzer mid- 
infrared detection. 

Based on these literature findings, we can estimate 
the expected number of planemos in NGC 1333. This 
is illustrated in Fig. [T2J In our survey area, we find 
19 brown dwa r fs, sey en of them previously identified by 
Wil king et all (|2004ft . In addition, the area covers 14 
more substellar objects from other sources (see Sect. l5.ip . 
With one exception, all these brown dwarfs exhibit opti- 
cal extinctions < 10 mag. For this extinction value, our 
survey is complete down to ~ 0.008 M@. The number 
of confirmed planetary-mass objects in the survey lim- 
its is zero. Based on the total sample covered in our 
spectroscopic observations, the la upper limit for this 
non-detection is about 2 objects. 

Assuming the planemo vs. brown dwarf ratio as seen in 
the a Ori cluster, we would expect to find 8 ± 3 planemos 
in NGC 1333, which has to be considered a lower limit, 
given the fact that the bottom of the mass function 
ha s not been found ye t in a Ori. The ONC survey 
bv lWeights et a l. (2009) predicts a number of planemos 
roughly comparable with the number of brown dwarfs, 
i.e. ~ 33 for NGC 1333. Based on the two other esti- 
mates, and ass uming a total cluster p opulation of ~ 150 
in NGC 1333 ([Gutermuth et al.l 12008ft we would expect 
< 20, most likely ~ 10 planemos in the cluster. Based 
on the mass function in Cha I, we expect 6-12 planemos 
in NGC 1333. Thus, all estimates for the expected num- 
ber of planetary-mass objects are as high as or higher 
than the 2a limit derived from our survey in NGC 1333. 
We conclude that NGC 1333 likely shows a deficit of 
planemos compared with all other star-forming regions 
where comparable survey depth has been reached. 

We arrive at essentially the same conclusion when we 
use temperatures instead of masses: The coolest object 
in our survey has an effective temperat ure of 2500 K, 
corre sponding to a spectral type of M9 (|Luhman et al.l 
2003). In crOri, th ere are 12 objects known with 
later spectral type s (jBarrado v Navascues et ah! 120011 : 
iMartm et all [2001): a significant fraction of the low- 
mass population is not spectros c opical ly confirmed. The 
ONC survey by I Weights et all (|2009ft gives 17 sources 
with >M9 and effective temperatures < 2500 K, about 
half of their total sample of substellar objects. In the 
Chal census, the four objects with planetary masses all 
have spectral types later than M9. In summary, work- 
ing with the maybe more robust effective temperatures 
instead of masses does not change the expected numbers 
of planetary-mass objects in any significant way. 

Our result i s sim ilar to one of the conclusions by 
iGreissl et all (|2007| ): Although they were sensitive 
to objects fainter t han those surveyed previously by 
Wilki ng" et al.l (|2004l ) , they did not detect lower mass ob- 
jects, only sources with higher extinction. We confirm 
this finding based on a much larger survey with a factor 



of ~ 300 more area coverage, encompassing the full ex- 
tent of the cluster. This result is even more remarkable 
when we take into account that the number of brown 
dwarfs in NGC 1333 is unusually high, as pointed out in 
Sect. El 

The best explanation for this outcome is a significant 
drop in the mass function in NGC 1333 around the Deu- 
terium burning limit, corresponding to effective tempera- 
tures below 2500 K in this cluster. Compared with other 
clusters like crOri and the ONC, the mass function of 
NGC 1333 does not extend down to planetary masses. 
As our survey is deeper than previous work in other re- 
gions and covers a substantial sample of objects with 
follow-up spectroscopy, this is the best evidence so far 
for a cut-off in the mass function at very low masses. It 
is plausible that our survey has detected the minimum 
mass limit for star formation in this particular cluster, 
at around 0.012-0.02 M Q . 

With the currently available sample, this is a 2a re- 
sult and thus needs further verification. If confirmed, it 
would point to strong environmental dependences in the 
mass function, particularly the mininum mass of star for- 
mation. In all other regions surveyed with comparable 
depth, the mass function extends well below 0.01 Mq, 
while the spectrum is truncated at ~ 0.012 — 0.02 Mq in 
our target cluster NGC 1333. 

At this point it is not clear what could cause a change 
in the minimum mass by a factor of two or more. Com- 
pared with regions as diverse as cOri, the ONC, and 
Chal, there is no evidence that NGC 1333 is differ- 
ent in any significant way in terms of initial conditions, 
dust properties, or metallicity. Based on the currently 
available census, the object density in NGC 1333 is 
slightly higher than in Cha I and significantly higher 
than in a Ori. While the number of low-mass mem- 
bers in all three regions is probably comparable and 
in the range of 150-250 objects, the radii of their dis- 
tribution are 0.3 pc in NGC 1333 , 0.5-0.7 pc in Cha I , 
and 3-5 pc in a Ori dSherrv et all l200l [Luhmanl [20071 
iGutermuth et al.l 12008ft . iBonnell et al.l (|2008ft recently 
suggested that brown dwarfs in clusters form through the 
gravitational fragmentation of infalling gas. This could 
explain the overabundance of substellar objects in the 
denser NGC 1333, but would also predict a higher num- 
ber of planetary-mass objects. On the other hand, the 
differences in density in those three regions are, as ex- 
plained in Sect. 14. 5[ most likely a consequence of pro- 
gressing dynamical evolution and thus a difference in age, 
and not necessarily in initial conditions. Thus, a consis- 
tent interpretation of our findings in terms of the forma- 
tion process is currently lacking. 

6. CONCLUSIONS 

We present the first results of the SONYC project, 
short for 'Substellar Objects in Nearby Young Clusters'. 
We have obtained deep optical and near-infrared obser- 
vations of the NGC 1333 cluster using the Subaru 8-m 
telescope. From the photometry in the i', z', J and K 
band, we have selected brown dwarf candidates. Follow- 
up infrared spectroscopy, also obtained at Subaru, is used 
to verify their nature. In the following, we summarize the 
most important findings of this survey: 

• In our large-scale, optical+near-infrared imaging 
survey, we reach completeness limits of 24.7 in i'- 
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Fig. 12. — Mass distribution of brown dwarfs in NGC 1333 (hatched histogram) and the deficit of planetary-mass objects. The three 
labelled datapoints show the predicted number of planemos in NGC 1333 based on the surveys in cOri, the ONC, and in Chal. Errorbars 
correspond to la. The x-values for these three datapoints roughly correspond to the completeness limits in the literature surveys for these 
regions. The non-detection of planetary-mass objects points to a deficit of objects in this mass regime compared with other clusters. The 
dotted lines refer to our detection limit for Av = 5 and 10 mag.. 



band, 20.8 in J-band, and 18.0 in K-band. In terms 
of object masses for members of NGC 1333, this 
corresponds roughly to mass limits of 0.008 M® 
for A v < 10 mag and 0.004 M Q for A v < 5 mag, 
based on the COND03 and DUSTY00 evolution- 
ary tracks. The spatial coverage of our survey 
is 0.25sqdeg and includes all previously identified 
cluster members as well as the central cloud cores 
NGC 1333-A and -B. 

• From the optical photometry, we select 196 ob- 
jects with i'-z' colors in the range as expected for 
substellar and planetary-mass cluster members, in- 
cluding a number of previously identified brown 
dwarfs. 36 of these objects are chosen for follow-up 
spectroscopy. This spectroscopic sample does not 
show any bias in spatial coverage or optical/near- 
infrared colors with respect to the full photometric 
candidate sample. 

• Using multi-object spectroscopy in H- and K-band, 
we confirm 28 objects close to the cluster center as 
late-type sources, 21 of them from the sample se- 
lected from photometry. This is based on water ab- 
sorption features in the H-band. The extinctions, 
derived from the J — K color, for the 28 objects 
cover the range from to 13 mag. Fitting model 
spectra to the observed ones, we find effective tem- 
peratures of 2500-3900 K with a typical uncertainty 
of ±200 K. 



• From the 28 confirmed late-type objects, 19 have 
effective temperatures of 3000 K or lower (and spec- 
tral types M6 or later). Combined with the clear 
indications for youth, this classifies them as proba- 
ble substellar members of NGC 1333. Of this sub- 
sample, seven have been previously identified by 
iWilking et all (pOOl . The remaining nine objects 
in the confirmed sample are possible steller mem- 
bers, which likely exhibit significant contamination 
by field dwarfs/giants. 

• The confirmed objects are strongly clustered 
around the core of NGC 1333. The majority 
of them shows evidence for disks in Spitzer mid- 
infrared colors. Both findings additionally provide 
evidence for membership. 

• Combining our survey results with previous stud- 
ies, the current census of spectroscopically con- 
firmed brown dwarfs in NGC 1333 is 33. For com- 
parison, the number of stellar members with masses 
between 0.08 and 1M Q is in the range of 50. This 
indicates a clear overabundance of brown dwarfs in 
NGC 1333; the ratio of substellar to stellar mem- 
bers with masses below 1 M Q is 1.5 ± 0.3, lower by 
a factor of 2-5 than in all other previously surveyed 
regions. 

• The low-mass limit of the confirmed brown dwarfs 
in NGC 1333 is 0.012-0.02 M©, but the complete- 
ness limits are at significantly lower masses. Thus, 



18 



Scholz et al. 



although being able to do so, we do not find a 
population of planetary-mass objects in this clus- 
ter. Scaling from literature results in other regions 
(a Ori, ONC, Cha I), we would expect to find 8- 
10 'planemos', but we find none. This indicates 
a cut-off in the mass spectrum around the Deu- 
terium burning limit in NGC 1333; we have possi- 
bly reached the minimum mass of the IMF. 

• If confirmed, these last two findings point to re- 
gional/environmental variations as a major factor 
in determining the number of very low-mass ob- 
jects and the minimum mass of the IMF. It is 
particularly difficult to explain simultaneously the 
overabundance of brown dwarfs and the deficit of 
planetary-mass objects in NGC 1333. Thus, our re- 
sults present a challenge for the current theoretical 
framework of star formation. 

• The mass spectrum in the brown dwarf and plane- 
tary mass regime necessarily relies on evolutionary 
models. Although we have made every effort to 
derive robust conclusions, we cannot get around 
the fact that evolutionary models lack reliable cal- 
ibrations at the lowest masses. This is particularly 
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true at young ages, where initial conditions may 
affect mass estimates, for example. Thus, the lack 
of objects with T < 2500 K in NGC 1333 could in 
principle be a sign of an unusual mass function or 
unusual initial conditions, or both. On the other 
hand, testing and calibrating the models require 
that we first identify objects at these masses and 
ages. Making relative comparisons between sub- 
stellar populations of different clusters could also 
help. Our SONYC project is an important step 
towards achieving these goals. 
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